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Abstract— In the case of transversely only loaded shallow shells, the nonlin- 
ear Donnell-Mushtari-Vlasov theory for large deflection of isotropic thin elastic 
shells leads to a system of two coupled nonlinear forth-order partial differential 
equations known as Marguerre's equations. This system involves two arbitrary 
elements - the curvature tensor of the shell middle-surface and the function of 
transversal load per unit surface area. In the present note, the point symme- 
try groups of Marguerre's equations are established, the corresponding group 
classification problem being solved. It is shown that Marguerre's equations are 
equivalent to the von Karman equations for large deflection of plates in the time- 
independent case and in the time-dependent case as well. It is also observed that 
the same holds true in respect of the field equations for anisotropic shallow shells. 



1 Introduction 

Within the framework of the nonlinear Donnell-Mushtari-Vlasov (DMV) theory (see, e.g., 
) the state of equilibrium of a transversely loaded thin isotropic elastic shell of uniform 
thickness is determined by the following system of two coupled nonlinear forth-order partial 
differential equations: 

DA 2 w - e a ^"w; a ^.,„ v ~ E a ^ v b a ^ v =p, m 
{l/Eh)A 2 <$> + (l/2)e a ^e^w., a/ 3W.^ + e a »e^b a pw^ v =0, U 

in two independent variables, associated with the coordinates on the shell middle-surface F, 
and two dependent variables - the transversal displacement function w, and Airy's stress 
function Here, e a P is the alternating tensor of F; b a p is the curvature tensor of F; D, 
E and h are the bending rigidity, Young's modulus and thickness of the shell, respectively 
(i.e., D, E and h are given constants); p is the function of transversal load per unit surface 
area; a semicolon is used for covariant differentiation with respect to the metric tensor a a /3 
of the surface F; A is the Laplace-Beltrami operator on F. Here and throughout, Greek 
indices have the range 1 , 2 and the usual summation convention over a repeated index (one 
subscript and one superscript) is employed. 

The present note is concerned with a special case of the nonlinear DMV theory; the 
so-called shallow shells are considered. In fact, this is an approximation of the theory, 
which is usually introduced as follows (cf., e.g. pi). Let ( X j X • X 3 ) be a fixed right-handed 
rectangular Cartesian coordinate system in the 3-dimensional Euclidean space in which the 
middle-surface F of a shell is embedded, and let this surface be given by the equation 

x 3 = f(x\x 2 ), {x 1 ^ 2 ) e n c R 2 , 

where / : R 2 — > R is assumed to be a single- valued and smooth function possessing as many 
derivatives as may be required on the domain fl. Let us take x 1 , x 2 to serve as coordinates on 
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the surface F. Then, relative to this coordinate system, the components of the fundamental 
tensors and the alternating tensor of F are given by the expressions: 

a a p = 6 afi +f, a f, P) b a p = ar 1 ' 2 f, aPi e a ^ a - 1/2 e^, (2) 

where 

a = det(a Q/3 ) = 1 + (/.i) 2 + (/. 2 ) 2 ; 

o~a/3 = <5 Q/3 is the Kronecker delta symbol; e"' 5 is the alternating symbol; here and in what 
follows, a comma is used for partial differentiation with respect to the coordinates on F. A 
shell is said to be shallow on the domain Qq C £1 when the inequalities 

\f, a \ \f,p\ < £ 2 < 1, s = const, 

hold for every point (x x ,x 2 ) E Qo- Hence, for shallow shells the quadratic terms in the 
right-hand sides of expressions (2) are small compared to unity and may be neglected. 
Thus, allowing for a relative error of order O (e 2 ) , one may regard the intrinsic geometry 
of the shell middle-surface F as Euclidean and (a; 1 , a; 2 ) may be thought of as an Euclidean 
coordinate system in which: 

O-a/3 = Saf3, (3) 
ba/3 = f,a/3, (4) 

e ap = e Q/3 ; (5) 

the mean curvature H of the surface F and its Gaussian curvature K (note that the latter 
is not necessarily equal to zero within the allowed relative error) take the form 

H = (1/2)5^1^, (6) 
K = (l/2)e a ^f ta pf <fiv , (7) 



and system (ED) reads 



D8 a ^w, aPlw - e^e^w^^ - e^e^b^^ = p, 
(l/Eh)6 a P6^®, afj ^ + {l/2)e a » e^w^pw^ + e^e^b af} w^ = 0, 



(8) 



Thus, we arrive at the equilibrium equations for shallow shells within the framework of the 
nonlinear DMV theory. Since equations (||) follow from Marguerre's shell theory || as well, 
they are also known as Marguerre 's equations for large deflection of plates with small initial 
curvature (i.e., shallow shells). These equations are well accepted and play an important 
role in the shell theory (see, e.g., [Q, || and the references therein). They also include as 
a special case, with b a p = 0, the well-known von Kdrmdn equations for large deflection of 
plates i]. 



2 Symmetry Groups 

The aim of the present work is to study, following , S and || , the invariance properties of 
system (||) relative to local one-parameter Lie groups of local point transformations acting 
on open subsets of the 4-dimensional Euclidean space R 4 , with coordinates (x 1 , x 2 , w, <£■), 
representing the involved independent and dependent variables. For that purpose Lie in- 
finitesimal technique is used and, as a rule, the results obtained are expressed in terms of the 
infinitesimal generators (operators) of the groups; in the present case, the latter are vector 
fields of the form 

where £ M ,?7 and ip are functions of the variables x , x , w and <&. The system considered 
involves an arbitrary tensor field - b a p , and an arbitrary function - p. This gives rise to a 
group classification problem with respect to the arbitrary element - the set {b a p,p}. 



The infinitesimal criterion of invariance leads to the following system of determining 
equations (DE system) for the components £ 1 ,£ 2 ,'r) and (p of the vector fields of form (^|) 
generating point symmetry groups admitted by system (^|) : 

? = C lX x + C 2 x 2 + C 3 , (10) 

£ 2 = -C2X 1 + Cix 2 + C 4 , (11) 

n = v (x\x 2 ), (12) 

tp = B 1 x 1 + B 2 x 2 + B 3 , (13) 

bant% + + ^Kp,n = -V,aP, (14) 

e^e^b^r),^ = 0, (15) 

DS°^5^ri, afilu/ = 2p^ + e"p, M , (16) 

where G\, C4, B\, B 2 , and B 3 are arbitrary real constants. This result is established in 
i jiol through the standard computational procedure (see J?], Sec. 5] or ||, Sec. 2.4]), and 
therefore we can assert that the system (Qj admits a vector field of form (||) if and only if 
- (16) hold. In other words, all possible symmetries of the kind under consideration 
inherent to system (ra) can be found via the solution of the DE system. So, the problem 
to solve consists in finding the solutions of this system. The main difficulty here comes out 
of the determining equations (]lj) - (|g|). They show that the space of solutions L to the 
DE system, i.e., the Lie algebra of symmetries associated with system (||), depends on the 
choice of the arbitrary element {b a f3,p}. At this juncture, we do face a group classification 
problem. That is to say, that we will have to determine all those specializations (special 
forms) of the arbitrary element for which system (^|) admits vector fields of the form (^J) . 

Proceeding to analyze this problem, we will first show that equations (|l2]), ( [l4| ) - ( |l6|) 
may be replaced by the following three equivalent ones: 

n = -e i U+A 1 x 1 +A 2 x 2 +A 3 , (17) 

where A\, A 2 and A 3 are arbitrary real constants, and 

2P^ + eP„ = 0, (18) 

2K^ + ^K tlt = 0, (19) 

where 

P = 2D6^H.^+p. (20) 

Indeed, ( |l7|) represents the general solution of equations ( |T^ ) of the form ( |l2|) when (fl ) , 
( [l0|) and (|ll[) hold; under the same assumption, by substituting ( |l7| ) into (15) and (16), 
and taking into account (0) and (0), after some algebra we obtain expressions (18) and (19). 
Using this result, hereafter we will assume that the DE system consists of equations (It), 
([OP, (0j, ®) and fifty. 

Now, a 6-dimensional space of solutions to the DE system arises immediately. It corre- 
sponds to the 6-dimensional Lie algebra Lq of the vector field whose components are given 
as follows 

e = 0, 

n = Aix 1 + A 2 x 2 + A 3 , (21) 
tp = Bix 1 + B 2 x 2 + B 3l 



where A%, A 2 , A3, B\, B 2 and B 3 are arbitrary real constants. Evidently, the above solu- 
tions of the DE system do not depend on the choice (specialization) of the arbitrary element. 



Hence, for any specialization of the arbitrary element system ^) admits the 6-parameter Lie 
group Gq generated by all linear combinations of the vector fields 



d 1 d 2 d d l d 2 d 

dw' x dw' x dw' a¥' x M' x a¥ ; 

these vector fields constitute a basis of the associated Lie algebra L . Simultaneously, Go is 
the largest Lie group of point transformations admitted by system |3) for any choice of the 
arbitrary element. Indeed, by setting 

/ = sin(x ) sin(cc 2 ), 

for instance, one can easily verify that for this particular choice of the arbitrary element the 
DE system has no other solutions besides those given by 

So far, we have obtained the so-called kernel of the full symmetry groups associated 
with system (||), that is the group Go- The next step is to identify the cases in which the 
system under consideration possesses larger groups of point symmetries. In the light of all 
the above, it means to characterize in a suitable manner all those specializations of the 
arbitrary element for which system (pf) admits vector fields of the form 



Y = ^— -eU-SZ (^0), (22) 



— -ei — 

dx^ dw 

where are given by the expressions ([!(]) and ([H]), keeping in mind that ( jTj| j and ( fil 
remain the only conditions (necessary and sufficient) for (0) to be invariant under a group 
generated by a vector field of form Note that all determining equations are thus taken 

into account. 

Taking into account , (Q) and ( pp| ) we can see at once that for / = p = the invariance 
conditions ( |l8| ) and ( pQj ) are satisfied. Hence, syste m ([§[ ), with / = p = 0, possesses a larger 
group of point symmetries (in addition to Go) and (|Io|) , (|Tl| ) and (|22|) show that this is the 
complete 4-parameter group of homothetic motions of the Euclidean plane. In this special 
case, coincide with the homogeneous von Karman equations, so that we have arrived 
at the result obtained in [0 (sec also p2| ). This example gives us a good motivation for 
studying the general case. 

We begin with the following observation. Given a vector field Y of form (p2j), the function 



w = w + f, (23) 

is an invariant of the corresponding Lie group of transformations acting on R 4 . Therefore, 
we can introduce new coordinates [x , x 2 , w, 3>) on R , the new dependent variable w being 
defined by (^3[), in which Y takes the form of an infinitesimal operator of the Lie algebra 
associated with the group of homothetic motions of the Euclidean plane, namely 

(note that the components i; 1 and £ 2 still have the form (|l7j) and (|ll]), respectively, as 
before), and system (||) reads 



DS a ^S^w, a ^ - e a »e^w, af3 <P^ = P, 
(l/Eh)6 a P6^$, af}fll/ + {l/2)e^e^w, a pw^ = K. 



(25) 



Thus, the problem of invariance of system (Sj with respect to a vector field (j^j converts into 
the problem of invariance of system ( jj^ ) under a vector field ( fj^| ) as a change of the variables 
does not affect the group properties of a system of differential equations. Now, taking into 
account the invariance conditions (|l8j) and ( |l9|) (note that they remain unchanged under the 
above coordinate transformation) , we can conclude that system (jjj[) admits a one-parameter 
group G of homothetic motions of the Euclidean plane generated by a vector field of form n24 



if and only if the corresponding arbitrary element {b a p,p} is such that the functions P and 
K, defined by formulae (Qj and ^idp, are invariants of G (when C\ — 0) or eigenfunctions 
(when C\ ^ 0) of its generator, the latter being regarded as an operator acting on the smooth 
functions £ : M — > R, M C R 2 . This result may be thought of as a general solution to 
the group classification problem under consideration in terms of the function p and the two 
characteristic invariants, H and K, of the shell middle-surface F. 



3 Equivalence Transformations 

Let us now briefly discuss, in the context of the shell theory, the meaning of the coordinate 
transformational : R 4 — » R 4 , [x , x 2 , w, $) i — ► (x , x , w + f, $) introduced in the previous 
Section. Omitting tilde's in equations (p5|), we can say that systems (||) and (pq ) belong 
to the same class - they have the same differential structure, and differ from one another 
only in the form of the arbitrary element. Hence, according to Definition 6.4], we may 
conclude that u> : R 4 — ► R 4 is an equivalence transformation for system^). In addition, 
we see that (^5j) is nothing but a system of nonhomogeneous von Karman equations with 
special right-hand sides. It is noteworthy that Marguerre 's equations (Q), i.e., the equilibrium 
equations for shallow shells, turned out to be equivalent to the von Karman equations, i.e., 
the equilibrium equations for plates, regardless of the invariance properties of system (^) . To 
the best of the author's knowledge, this fact has not been noticed before in the literature, 
though both the von Karman equations and Marguerre's equations have been studied and 
utilized for many years until now (see, e.g., 0, g], Q, ||, pT| , |12 , jl3| and the references 



therein). In our opinion, the established correspondence between these two systems of 
equations will certainly be of use in solving a wide range of problems arising in the shell 
theory. In particular, the results for the von Karman equations obtained in Q can be easily 
conveyed to the theory of shallow shells. 

Several additional applications of the transformation to can by given too. Evidently, it 
maps: 

• the time-dependent Marguerre's equations 

DS a ^S^w, a ^ v - e^e^w,^,^ - e a ^e^b a ^,^ + pw = p, 

q/3m „ + (l/2)e a " e^w^w^ + e a »e^b a(i w^ v = 0, 

(where p is the mass per unit area of the shell middle-surface, and a superposed dot is used 
to denote partial derivative with respect to the time t) into the time-dependent von Karman 
equations 

D8 a n^w, a ^ v - e^e^w,^,^ + pw = P, 
{l/Eh)8 a H^^, a0ll u + (l/2)e a »e^w, af 3W^ = K; 

• the time-dependent equations for anisotropic shallow shells 

D a ^w. a[3 ^ - e^e^w^^ - e^e^b^^ + pw = p, 
E a ^^, a0 ^ + (l/2)e a »e^w, a pw^ + e^e^b^w.^ = 0, 

(where D a ^^ v and E -^^ denote the material constants) into the following system of von 
Karman-type equations 

D af3 ^w^ - e a ^> a/ 3$ ^ + pw = D a ^f. afj ^+p, 
£ Q/3 ^$ + {l/2)e^e^w^w.^ = K. 

The above list could be continued. 
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